GONG p-mode frequency changes with solar activity 
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ABSTRACT 

We present a correlation analysis of GONG p-mode frequencies with nine 
solar activity indices for the period 1995 August to 1997 August. This study 
includes spherical harmonic degrees in the range 2 to 150 and the frequency 
range of 1500-3500 /iHz. Using three statistical tests, the measured mean 
frequency shifts show strong to good correlation with activity indices. A 
decrease of 0.06 /iHz in frequency, during the descending phase of solar cycle 22 
and an increase of 0.04 /zHz in the ascending phase of solar cycle 23 is observed. 
These results provide the first evidence for change in p-mode frequencies around 
the declining phase of cycle 22 and beginning of new cycle 23. This analysis 
further confirms that the temporal behavior of the solar frequency shifts closely 
follow the phase of the solar activity cycle. 

Subject headings: Sun: activity Sun: oscillation 
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1. INTRODUCTION 

Solar surface phenomena such as sunspots, magnetic field, 11-year solar cycle, and 
associated activities are the consequence of dynamical processes occurring inside the Sun. 
The fundamental parameters responsible for most of the activities are the magnetic and 
velocity fields, thus the acoustic mode (p-mode) frequencies could be influenced by the 
solar activity. Hence, a relation between the solar frequencies and activity indices may 
throw some light on the changes occurring deep in the solar interior and provide clues 
to the mechanism of solar cycle. The first evidence that the frequencies of the p-modes 
change with the solar cycle, came from the observations of low degree (£) modes. From the 
low-£ ACRIM data a decrease in the frequencies of 0.42 /xHz between 1980 and 1984 was 
reported by Woodard and Noyes (1985). It was suggested by |Kuhn ( 1 988 )| that the observed 



variation in frequencies is due to the changes in the internal temperature structure and are 
correlated with the variation in the surface temperature. Subsequently, |Goldreich et al. 



1991)| theoretically showed that the variations in the solar p-mode eigen frequencies are 



related to the perturbations in the magnetic flux at the suns surface. Observations made 
during the rising phase of the solar cycle 22 by the Big Bear Solar Observatory (BBSO) 
group found an increase of about 0.4 /zHz in mean frequency during the period 1986-89 
( [Libbrecht fc Woodard 1990| ). Further, it was shown that the frequencies varied on a time 
scale as short as 3 weeks with the absolute value of the magnetic field strength over the 
visible solar disk (Woodard et al. 1991). These results confirmed the theoretical work by 
Goldreich et al. ( 1991 )| . Now, there is growing evidence that the mode frequencies vary with 



solar activity. Bachmann fc Brown (199*3] using the High Altitude Observatory's (HAO) 



Fourier Tachometer data for I between 20 and 60 in the frequency range of 2600-3200 /iHz 
showed that the p-mode frequency shifts correlate remarkably well with six solar activity 
indices during 1984 October through 1990 November. For the low degree modes (£ < 3) 
and during different time epochs, several groups ( [Anguera Gubau et al. 1992] , Regulo et al. 
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1994, Jimenez-Reyes et al. 1998 ) have shown that the frequency shifts are well correlated 



with solar activity cycle. 

It has also been shown that the frequency splitting coefficients vary with solar cycle 
(Kuhn 1988, [Woodard and Libbrecht 1993| ). It was pointed out that the even-order 



coefficients are sensitive to the solar cycle or to the latitude-dependent properties while 
the odd-order coefficients reflects the advective, latitudinally symmetric part of the 
perturbations caused by rotation. Recently, using the GONG data, |Anderson, Howe,~^ 



[Komm (1998)| have found evidence for small shifts in the central frequencies and splitting 



coefficients. From the solar oscillation data obtained from Michelson Doppler Imager on 
board SOHO, during 1996 May 1 through 1997 April 25, Uzicmbowski et al. (1998) also 
detected a significant trend in the splitting coefficients. 

Motivation for the present work arose due to the availability of high precision 
frequencies from the GONG network, and for extending the earlier studies to include the 
declining phase of cycle 22 and the beginning of cycle 23. In this study, we also look for a 
possible correlation between the p- mode frequency shifts and nine solar activity indices 
representing the photospheric, chromospheric and coronal activities. 



OSCILLATION DATA AND SOLAR ACTIVITY INDICES 



The GONG data ( [Hill et al. 1996| ) used in this study consist of 8 data sets covering a 



period of two years from 1995 August 23 to 1997 August 11, in the frequency range between 
1500 to 3500 //Hz and i from 2 to 150. This period covers the GONG months (GM) 4 to 9 
(1995 August 23 to 1996 March 25) each of 36 days, GM 12-14 (1996 June 6-September 21) 
and GM 21-23 (1997 April 26-August 11) of 108 days each. The period from 1995 August 
through 1996 September covers the declining phase of the solar cycle 22, whereas the period 
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from 1997 April-August refers to the beginning of cycle 23. 

For calculating the frequency shifts, we have used the frequency of GM 12-14 as a 
reference. The choice of this standard frequency has been made in view of two reasons; 
(a) that this datum lies in the middle of the period covered in this study, and (b) that it 
represents the period of the minimum solar activity. The mean frequency shift is found 
from the relation: 

M*) = £%W)/£% (i) 

where Q n £ is the inertia ratio as defined by |Chistensen-Dalsgaard fc Berthomieu (199T) . 
a n £ is the error in frequency measurement, as provided by GONG and 8vnz(\) is the change 
in the measured frequencies for a given £, and radial order, n. We have neglected those 
mode frequencies for which a n £ > 0.1 /iHz, to avoid large errors in calculating the mean 
shifts. The resulting mean weighted frequency shift 8v is plotted against the GONG months 
for four different £ ranges (2 < £ < 19, 20 < £ < 60, 61 < £ < 150, and 2 < £ < 150) 
and is shown in Figure 1. For the interval GM 4 to 9, the frequency shift shows a systematic 
decrease of about 0.06 /zHz for all £ values while for GM 9 to 23, the mean shift shows 
an increase of 0.04 /zHz, indicating that the solar oscillation frequency varies with solar 
activity. 



EDITOR: PLACE FIGURE | HERE. 



We have correlated the mean frequency shifts with different solar activity indices: the 
International sunspot number, Rj, obtained from the Solar Geophysical Data (SGD); 
KPMI, Kitt Peak Magnetic Index from Kitt peak full disk magnetograms ( [Harvey 1984|) ; 
SMMF, Stanford Mean Magnetic Field from SGD; MPSI, Magnetic Plage Strength Index 
from Mount Wilson magnetograms ( Ulrich et al. 199~I| ); FI, total flare index from SGD and 
Atag 1999]; He I, equivalent width of Helium I 10830A line, averaged over the whole disk 
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from Kitt peak; Mg II, core-to- wing ratio of Mg II line at 2800 A from SUVSIM (SGD); F 10 
integrated radio flux at 10.7 cm from SGD; Fe XIV, Coronal line intensity at 5303A from 
SGD. A mean value was computed for each activity index over the interval corresponding 
to the same GONG month. 



3. ANALYSIS AND RESULTS 

To study the relative variation in frequency shift 5v with activity index i, we assume a 
linear relationship of the form: 

8u = ai + b, (2) 

where 5v includes the mean error in the measured frequencies. The slope, a and the 
intercept, h are obtained by performing a linear least square fit. Figure 2 shows a typical 
plot of our analysis for all modes between 2 to 150. The solid line represents the best 
fit regression line and shows that the data is consistent with the assumption of a linear 
relationship. The bars represent la error in fitting. We also carried out x 2 -test which takes 
into account the statistical uncertainties of each of the frequency shift. As there is no 
reliable estimate available of the uncertainties in the measurement of activity indices, those 
are not included in the fitting. 

EDITOR: PLACE FIGURE g HERE. 



EDITOR: PLACE TABLE HERE. 



The linear relationship given in equation (2) is further tested by calculating the 
parametric Pearson's coefficient, r p , and the non-parametric Spearman's rank correlation 
coefficient, r s along with their two-sided significance, P p , and P s respectively, for all the 
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£ ranges. As a typical example, these parameters are given in Table 1, for the £ range of 
2 to 150. From the table, it is clear that a positive correlation exists for all the activity 
indices. We note that the best correlation is obtained for Mg II index, while Fio and 
Fe XIV indices show good correlation. SMMF show poor correlation, perhaps due to large 
gaps in the available data. Normally, we expect better correlation should exist between the 
surface magnetic activity and Su, as this is the fundamental parameter for all solar activity. 
However, from Table 1, we notice that the radiative indices (-Flo, Mg II and Fe XIV) show 
better correlation as compared to magnetic field indices. Backmann and Brown (1993) 
using the HAO data had obtained a similar result. The reason for such a difference in the 
level of correlation needs further investigation. 

In order to investigate the degree dependence of frequency shifts with activity indices, 
we carried out a statistical analysis for three different ranges of £. It is seen that the fitting 
parameters, a and b do not show significant variation with t. The lowest \ 2 values are found 
for 2 < £ < 19, while for 20 < £ < 60 and 61 < £ < 150 the x 2 values are higher. It is 
further noted that the correlation coefficients, r p , and r s for all £ ranges have values between 
0.75 to 1.0 for all solar indices except for He I and SMMF. Comparing r p and r s values for 
all activity indices, it is noted that for 2 < £ < 19 these values are systematically less, 
indicating a week correlation as compared to the other ^-ranges. Thus, there appears to be 
some evidence that the frequency shift depends on £. 

EDITOR: PLACE TABLE | HERE. 



Comparing our results for 20 < £ < 60 with those of |Woodard et al. (1991)1 , and 



Bachmann fe Brown ( 1993 )| for modes of similar degrees (Table 2), we find that the 



magnitude of slope a, is comparable in all the three cases for the KPMI index, the only 
index common in all the three analysis. However, a small difference of 0.003 /iHz between 
the GONG and BBSO data and 0.001 fiRz between the GONG and HAO data is noticed. 



- 8- 



This difference may be due to different frequency intervals and phase of the solar cycle. A 
detailed comparison between BBSO and HAO results was made by [Regulo et al. (1994) . 
The other four indices considered in HAO and our work show similar correlation. 

EDITOR: PLACE FIGURE g HERE. 

The temporal behavior of the solar frequency shifts and the mean monthly sunspot 
number R z is shown in Figure 3. The diamonds represent the HAO data for the period 
1984 October to 1990 October, while the triangles indicate the GONG data for the period 
1995 August to 1997 August. A polynomial of order n of the form 

n 

5u = ao + a id % (3) 

i=l 

is fitted separately to each set of frequency shift 8v since the reference frequencies for HAO 
and GONG data are different. Here, a®, and a* are the fitting coefficients and d the number of 
days after 1981 January 1. The best fit for HAO data is obtained with a polynomial of order 
4, while a polynomial of order 2 fits the GONG data. The fitting coefficients for HAO data 
are; a = 0.557, ai = 1.583xl0~ 4 , a 2 = -1.107xl0" 6 , a 3 = 5.839xl0- 10 , a 4 = -8.282xl0~ 14 , 
and for GONG data are; a = 10.948, a x = -3.847xl0" 3 , a 2 = 3.369xl0~ 7 . From this plot 
it is evident that the mean frequency shifts systematically follow the solar activity cycle 22 
and the beginning of cycle 23. 

We conclude that a positive and linear correlation exists between the frequency shift 
and the nine solar activity indices. It is observed that the radiative indices show better 
correlation as compared to the magnetic indices. This analysis shows a decrease of 0.06 /xHz 
in weighted frequency during the declining phase of solar cycle 22 and an increase of 0.04 
/iHz during the ascending phase of cycle 23. Our finding confirms that the temporal 
behavior of the solar frequency shifts, closely follow the phase of the solar activity cycle and 
we conjecture that the same trend will continue during cycle 23. 
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Fig. 1. — Mean frequency shifts for p-modes weighted by the mode inertia in the frequency 
range of 1500-3500 //Hz for different ranges of £; 2 < £ < 19 (dotted line), 20 < £ < 60 
(dashed), 61 < £ < 150 (dash-dot), and 2 < £ < 150 (solid). Error bars represent mean 
error in shift. 

Fig. 2. — Comparison of weighted mean frequency shift with nine different activity indices 
in the £ range of 2 to 150. These are (a) sunspot number, (b) Kitt Peak magnetic index 
(Gauss), (c) 10.7 cm flux (sfu), (d) He I 10830 A equivalent width (mA), (e) Stanford mean 
magnetic field (/iTesla), (f) Magnetic plage strength index (Gauss), (g) Total flare index, (h) 
Coronal line intensity (10 16 W/sr), (i) Mg II core-to- wing ratio. The solid line represents the 
linear fit. The error bars indicate la error of the fitting. 

Fig. 3. — The top panel shows the time variation of frequency shifts for the HAO (diamond) 
and GONG (triangles) data for a total period from 1984 to 1997. The best fit for GONG 
data is obtained with a polynomial of order 2 while the HAO data is best represented by a 
polynomial of order 4. For comparison, the bottom panel shows the time variation of the 
mean monthly sunspot number R z . Error bars for GONG data are not clearly visible as 
these are very small. 
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Table 1. Fitting and correlation statistics for 2 < i < 150 



Activity index 


slope a 


intercept b (/iHz) 


x 2 






r s 


Ps 


Ri 


0.002 ± 0.0001 /iHz 


-0.034 ± 0.001 


344 


0.82 


0.023 


0.93 


0.002 


KPMI 


0.029 ± 0.001 /iHz/G 


-0.221 ± 0.008 


461 


0.81 


0.025 


0.71 


0.071 


Fio 


0.005 ± 0.0002 /iHz/sfu 


-0.345 ± 0.011 


144 


0.93 


0.002 


0.93 


0.002 


He I 


0.010 ± 0.0004 /iHz/mA 


-0.449 ± 0.017 


480 


0.76 


0.049 


0.61 


0.148 


SMMF 


0.0004 ± 0.0002 fiRz/f/T 


-0.009 ± 0.002 


1130 


0.16 


0.729 


0.21 


0.645 


MPSI 


0.191 ± 0.006 fiRz/G 


-0.031 ± 0.0009 


266 


0.87 


0.011 


0.93 


0.002 


FI 


0.023 ± 0.0008 /iHz 


-0.017 ± 0.0005 


373 


0.83 


0.020 


0.93 


0.002 


Fe XIV 


0.012 ± 0.0004* 


-0.048 ± 0.001 


128 


0.93 


0.006 


0.94 


0.005 


Mg II 


38.210 ± 1.172 /iHz 


-9.712 ± 0.297 


68 


0.97 


0.003 


0.90 


0.005 



>Hz sr/10 16 W 
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Table 2. Comparison between GONG and HAO results a 



GONG analysis HAO analysis 



Activity Index 


a 


b 


x 2 


a 


b 


x 2 


KPMI* 


0.024 


-0.182 


216 


0.023 


-0.306 


434 


Fio 


0.004 


-0.307 


60 


0.003 


-0.345 


153 


He I 


0.009 


-0.410 


170 


0.011 


-0.667 


184 


MPSI 


0.171 


-0.031 


94 


0.136 


-0.168 


166 


Mg II 


33.986 


-8.643 


22 


20.96 


-5.700 


134 



a Units of a and b are same as given in Table 1 
*for BBSO, a = 0.027 /iHz/G 
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